Single photons from coupled quantum modes 
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Single photon emitters often rely on a strong nonlinearity to make the behaviour of a quantum 
mode susceptible to a change in the number of quanta between one and two. In most systems the 
strength of nonlinearity is weak, such that changes at the single quantum level have little effect. 
Here, we consider coupled quantum modes and find that they can be strongly sensitive at the single 
quantum level, even if nonlinear interactions are modest. As examples, we consider solid-state 
implementations based on the tunneling of polaritons between quantum boxes or their parametric 
modes in a microcavity. We find that these systems can act as promising single photon emitters. 



Introduction. — The construction of single photon 
sources P, 0] is a current aim of quantum nonlinear op- 
tics. Aside contributing to the security of quantum cryp- 
tography single photon sources are useful elsewhere, 
for example in schemes for quantum computation using 
only linear optics and photodetection For some ap- 
plications it is enough to reduce the intensity of a laser 
source to obtain single photons with a probability limited 
by the Poisson distribution. To do better than Poisson 
statistics one requires some form of nonlinearity. How- 
ever, when one works in the single photon regime a strong 
nonlinearity is not so easy to find. 

In semiconductor microcavitics, light is strongly cou- 
pled to quantum well excitons resulting in new quasipar- 
ticlcs known as polaritons. Taking the best from both 
parents, polaritons have attracted particular attention 
for over a decade due to their strong nonlinearity (inher- 
ited from excitons) as well as their fast dynamics, long co- 
herence and ability to couple to external light (features of 
photons). Polariton-polariton interactions have resulted 
in micron-sized optical parametric oscillators Q , opti- 
cal gates Q , spontaneous coherence , low threshold 
lasing at room temperature 12 - 3] and superfluidity fl5| . 
Whilst these effects involve many polaritons at once, we 
wish to focus on the single quantum regime. In planar 
cavities, quantum effects such as squeezing have been re- 
ported and several studies on quantum correlations un- 
dertaken 1(1- III- More pronounced effects at the single 



polariton level are expected in quantum boxes [20H22I ]. 
where polaritons are fully confined in three-dimensions 
and forced to interact even more strongly. Available re- 
cently, such confinement has encouraging prospects for 
single photon sources. 

It has been predicted that for a very strong nonlinear- 
ity, the presence of a single polariton can block the res- 
onant injection of another [23j . analogous to the photon 
blockade 24| of nonlinear cavities. However, to obtain 
a strong enough nonlinearity for a single photon source, 
an extremely small quantum box is required (with size of 
the order of 200nm). Although one may anticipate such 
a system in the future, current systems do not display 
such a strong nonlinearity - whilst high nonlinearity is 



present in semiconductor microcavities, the energy shift 
caused by two interacting polaritons remains small. 

We consider theoretically two coupled quantum boxes 
and show that the coupling can dramatically enhance the 
characteristics of single photon devices. By solving the 
quantum master equation for the density matrix, we find 
strong single photon statistics for values of the polariton- 
polariton interaction strength corresponding to today's 
systems. We show that this is due to correlations be- 
tween the quantum fluctuations in the two boxes, allow- 
ing a much stronger sensitivity of the system to the pop- 
ulation compared to the single mode case. We expect 
that the coupling under study can also be exploited in 
analogous systems such as coupled nonlinear cavities or 
coupled photonic crystal cavities (25|. Finally, we show 
how mode coupling in parametric oscillators can also en- 
hance single photon statistics due to selection rules. 

A pair of linearly coupled modes. — Consider a pair of 
quantum modes described by creation operators a\ and 
a\ respectively. As an example, we imagine the lowest 
energy polariton modes of two spatially separated micro- 
cavity quantum boxes. In each box, polariton-polariton 
interactions are characterized by an interaction strength 
a. The boxes are spatially separated such that there 
are no significant nonlinear interactions between boxes. 
However, the boxes are close enough together such that 
particles can tunnel from one box to the other, at a rate 
given by the tunneling constant J. The Hamiltonian is: 
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where E\ and E 2 are the uncoupled energy levels of the 
two quantum modes and F represents a coherent excita- 
tion of the first mode. With quantum boxes, this would 
be a laser excitation focused onto the first quantum box. 
We define our energy scale such that the pump energy is 
zero. The evolution equation of the corresponding den- 
sity matrix, p, is: 
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where the last term represents the standard Lindblad dis- 
sipation characterized by decay rate T. Equation ([2]) can 
be solved by expanding the density matrix over a particle 
number basis in a similar way to that done in Ref. [23|; 
one truncates at a given particle number and propagates 
in time from the vacuum to the steady state [26[. 

For a pair of quantum boxes of 3/im size, separated 
by IjUto, a typical value of the tunnel constant is J = 
0.5meV. Since J > Y, strong coupling takes place and 
the single particle eigenmodes are the symmetric and an- 
tisymmetric modes spanning the two wells 1271. We take 
a = 0.012meV, a value measured in Ref. [28[ for con- 
densed polaritons occupying spot sizes of ~ 3fim. Al- 
though the pump acts directly on the first well it effec- 
tively pumps the second well due to an interference effect. 
This can be understood by considering the Heiscnberg 
equations for the symmetric (a+ = sin </>&i +cos 4>d 2 ) and 
antisymmetric (a_ = cos^ai — sin0a2) field operators, 
with cigenencrgies E+ and E— respectively: 



E + a+ + -Fshw 



da— 
~aT 



E-d- + Fcos< 



For E\ « Ei, (j> « 7r/4 and in the steady state one finds 
that a± = Since E+ and E- have different signs, 

d + and a_ are excited with different signs such that a 2 
is excited instead of a\. 

A key quantity in quantum optics is the second order 
correlation function, defined as: 
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When n and m correspond to the same mode and t = t', 
this quantity also measures the performance of a single 
photon source; an ideal source has 52,7m (0) = 0, whilst a 
classical source has 32,nn(0) = 1- 

For fixed pump intensity, the dependence of 32,11(0) 
on the energy levels of the two wells is shown in Fig. 
The optimum (smallest) (72,11(0) was attained when E\ = 
0.07meV and E 2 = 0.05mcV. The 32,11(0) depends 
mostly on the energy of the second well, E2 , and this vari- 
ation is shown again in Fig. [T}) along with the variation 
°f 32,22(0), and the average well populations, (Ni) and 
(JVjj) respectively. Unlike in Fig. [TJi, the pump energy is 
also varied to maintain a constant detuning between the 
pump and the lowest energy (symmetric) single particle 
eigenstate. This allows a better test of the variation of E2 
since the average populations do not change drastically 
and reveals that whilst the second well has 32,22(0) ~ f , 
varying E 2 has a dramatic effect on 32,11(0). 

Utilizing the first well as a single photon source, one 
finds that at a pump amplitude giving (N%) = 0.02 the 
probability of having more than one photon is 0.18%; 
this is five times better than the failure rate of devices 
based on spontaneous parametric down conversion 
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FIG. 1: a) Variation of ff 2 ,n(0) with E 1 and E 2 . b) De- 
pendence of the equal time correlation function and average 
populations on E 2 for E 1 — 0.07meV. c) 32,11 (t) for the op- 
timum parameters from (a), d) (72,11(0) from Eq.((5| (solid) 
and comparison to the single mode case (dotted). In all panels 
T = 0.2meV and in (a-c) F = O.lmeV. 



Spectral filtering could further reduce emission from the 
m > 2 states, providing extra improvement. Fig. [TJ: 
shows the unequal time second order correlation func- 
tions [2^] , which oscillate at half the Rabi oscillation pe- 
riod arising from the J = 0.5meV coupling. 

To better understand the low value of 32,11(0) we car- 
ried out analytical calculations, extending the method of 
Ref. [2^], which applies directly to the single mode case. 
We use stochastic (Langevin-type) equations for the evo- 
lution of quantum fields (30j . Lowest-order fluctuations 
of the fields around their mean values can be found by 
solving the linearized (a n i-> d n + 5a n ) version of the 
equations (26[ |. Choosing the convention that d\ is real, 
the 2 nd order correlation can be written 130| : 



32,11 = 1 + — [(SalSai) + {(6aiS ai )}] (4) 
ni 



which yields [26( : 
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where ( = Y/ (2(E\ + 4ani)). The mean field values, 
01, can be obtained as in Ref. 27j. To calculate all sec- 
ond order correlations we have extended the method of 
Ref. [3l|. In Fig. QJi we compare the result of Eq. [S] 
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for the two mode (solid line) and single mode (dotted 
line) cases. The single mode value is obtained by setting 
J = and matches the result from Ref . [29[ . For parame- 
ters corresponding to the optimum g2.1i! the last term in 
Eq. [5] makes a strong negative contribution [26|. In other 
words, the correlated noise fluctuations (60,160.2) drive 
the low value of 52,11 • These correlations are a result of 
the interplay between nonlinearity (in the limit a \-> 0, 
(601602) 1— ► 0) and tunnelling. This is a very different 
mechanism from the polariton blockade j29j , which orig- 
inate from the third term in Eq. [5] that vanishes in the 
present regime where ani t 2 <C T. 

For pioneering experiments one may also consider re- 
placing the coupled quantum box modes by the circu- 
larly polarized spin modes of a single quantum box or a 
localized state in a planar microcavity. Magnetic fields 
parallel and perpendicular to the growth direction would 
allow the tuning of the spin energy levels (as suggested as 
a control method of the original polariton blockade [32I ]) 
and the coupling constant J, respectively. 

Parametrically coupled modes. — Given the attention 
devoted to parametric processes in microcavities at the 
beginning of the millennium , it is only natural for 
us to ask whether sub-Poisson statistics can also be de- 
rived from pair scattering processes. Again, a variety 
of systems can be imagined, including the eigenmodes 
of a localized quantum box |33[. The Hamiltonian now 
includes three separate modes characterized by creation 
operators a\, a 2 and a 3 respectively: 

H = ^2 \ End^On + aal l al l a, n a n + 2a al l al n a n a m 

n—l \ m^n J 

+ 2a faf ajc^c^ + a\a\aia^j + Fa\ + F*ci2 (6) 

The scattering terms represent the same selection rules 
as when one deals with a pump mode that can scatter in 
pairs to signal and idler modes, as in intra-branch [D, Q 
and inter-branch Q scattering in planar microcavities. 
The Hamiltonian can be diagonalized on a number state 
manifold in the absence of the pump terms. The en- 
ergy levels are shown in Fig. [2j The only state in the 
ni = 2 manifold that can pollute the value of 52,11 (0) is 
the lowest lying 1 200) state. An analysis of the allowed 
transitions shows that this state can only be reached from 
decay of an ni — 3 state. In the low occupation limit the 
system is not expected to visit the n\ = 3 manifold very 
frequently, leading us to expect a low value of 52,11(0). 

Using a similar evolution equation to Eq.Q, 52,11(0) is 
calculated in Fig. [3^, for different values of the energy lev- 
els, E2 and £3. A study varying also Ei, showed that the 
optimum parameters are for the case when the pumped 
mode is at an energy slightly higher than the pump en- 
ergy, E 2 = 0.12meV and the sum of the signal and idler 
mode energies is resonant with the pump (see Fig. [SJd) . 
In fact 52.11 (0) is found only to vary if Ei +E3 is changed, 
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FIG. 2: Eigenstates from the diagonalization of the Hamilto- 
nian Eq.® without the pump terms on the particle num- 
ber manifold. Parameters: E\ = OmeV, E\ = 0.12meV, 
E 2 = OmeV, a = 0.012meV, Y = O.lmeV. The transitions 
caused by optical pumping and decay are shown by the solid 
and dashed arrows respectively. 
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FIG. 3: a) Variation of 32,11(0) with E 2 and E 3 for E\ = 0, 
F — O.lmeV. b) The optimum energy levels, c) 32,11 for 
the optimum parameters. Other parameters: T = O.lmeV, 
a — 0.012meV. d) Comparison of the 32,11(0) available from 
the single mode polariton blockade (with pump resonant with 
the bare mode energy), the two coupled mode case (from 
Fig[T|) and the three parametrically coupled modes. 32,11(0) 
is plotted as a function of the average occupation number of 
the signal mode. Dashed curves show the effect of dephasing. 
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that is, a selection of values of E\ and can be used to 
obtain optimum results. 

The time dependent second order correlation is shown 
in Fig. [3b. For zero delay we obtain 52,11(0) = 0.28. It is 
important to remember that we are working with a value 
of the polariton-polariton interaction strength available 
in two dimensional planar cavities (28[; much lower val- 
ues of 32,11(0) would appear in more confined systems in 
which the strength of interactions is higher. 

For comparison, the second order coherence function 
is shown in Fig. [3Ji for three cases using the same value 
of the interaction strength: the single mode polariton 
blockade [IH; the two coupled well case (studied in the 
previous section); and the case of three paramctrically 
coupled modes. It is clear that the use of schemes in- 
volving two or three coupled modes can give statistics 
closer to that of a single photon device than the single 
mode case. Whilst in all schemes the value of 32,11(0) de- 
creases with the average population of the signal mode, 
for the single mode and three mode case 32,11(0) tends 
to a constant value as (N) is decreased. 

Finally, we have considered the effect of pure dephasing 
(associated with exciton-phonon scattering) by adding 
the term: 
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to the right-hand side of Eq.©. Taking T P = 0.3/ieV 
(an upper estimate of the dephasing from Ref. (HI) gives 
the dashed curves shown in Fig. (3Ji. 

Conclusion. — We considered the use of a pair of cou- 
pled quantum boxes as a single photon source. We 
present the general idea that coupling can dramatically 
improve the single photon statistics compared to the 
single-mode case, through noise correlations. With com- 
petitive characteristics one may choose quantum boxes in 
a solid-state system, which offer fast (picosecond scale) 
relaxation rates; compact size; emission into a well- 
defined spatial mode (a luxury not always present when 
working with quantum dots); and wavelengths compat- 
ible with transmission through silica fibres and photo- 
detection with silicon based photon-counters. Alterna- 
tively, one may consider coupling nonlinear cavities or 
using paramctrically coupled modes in planar microcav- 
ities. Indeed we anticipate more studies specific to each 
system and experimental verifications. The nonlinear- 
ity is available with our present technology and we have 
several paths to choose from. 

We thank M. Wouters and A. Fiorc for useful discus- 
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